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Level-spacing distributions of the Gaussian Unitary Ensemble (GUE) of random matrix theory are expressed 
in terms of solutions of coupled differential equations. Series solutions up to order 50 in the level spacing 
are obtained, thus providing a very good description of the small-spacing part of the level-spacing distri- 
bution, which can be used to make comparisons with experimental or numerical data. The level-spacing 
distributions can be obtained by solving the system of differential equations numerically. 
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1 Introduction 

Ever since the pioneering works of Wigner and Dyson (l]0|3]|4l|5l|6), random matrix theory has been a 
major tool in the investigation of complex systems in physics, see e.g. Refs. QUI- In particular, the three 
"classical" random matrix ensembles, the Gaussian Orthogonal Ensemble (GOE), the Gaussian Unitary 
Ensemble (GUE) and the Gaussian Sympletic Ensemble (GSE) have been studied in much detail 1 9 1 and 
have proven of relevance in various applications, not alone in physics, see e.g. | 1 1 for examples. 

One of the characteristic features that is frequently used in investigations is the distribution of differ- 
ences between eigenvalues. These level-spacing distributions are universal quantities, in the sense that 
they correctly describe the level-spacing distributions in many complex systems. In many applications, 
experimentally observed or numerically calculated discrete data, such as energy levels of large nuclei, are 
analyzed and their spacing distribution is compared to that of the appropriate random matrix ensemble, 
where the choice of the appropriate ensemble is dictated by the symmetry of the system. Besides the spac- 
ing between neighboring levels, also the spacing between levels separated by a fixed number of levels, 
i.e., next-nearest neighbors, next-next-nearest neighbors, etc., has been considered. A recent example of 
an astonishingly precise agreement between energy levels of a Hamiltonian and random matrix distribu- 
tions was found in a simple tight-binding model of an electron moving on a planar quasiperiodic graph 
II 11 1 1 21 1X31 1 1 41 . which is used to model electrons in quasicrystals 1 151 1161 . But the distribution functions 
of random matrix theory are even of relevance in pure mathematics; the paradigm of an application in 
number theory is the distribution of zeros of Riemann's zeta function on the critical line 1171 . 

However, there exists no simple closed form for the level-spacing distributions of the above-mentioned 
ensembles of random matrix theory. Thus, in many cases, the empirical data are actually compared to a 
so-called Wigner surmise, which corresponds to the spacing distribution of neighboring levels obtained 
from two-by-two random matrices, rather than from the limiting distribution for matrices of infinite size. 
For example, in the GOE case, the difference between the Wigner surmise and the true GOE spacing 
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distribution is rather small, but the example of the quasiperiodic tight-binding model II II 1121 1731 1141 
shows that it may indeed be detected in physical applications. For the GUE case, the situation is somewhat 
better as the complete small spacing expansion for the nearest neighbor spacing distribution is known |7 1. 
Numerical values, series expansions and a Pade analysis for the nearest-neighbor spacing distributions can 
be found in 191 171 1181 . However, level spacings beyond the nearest-neighbor distributions, the spacing 
distributions of levels separated by a number of other levels, have scarcely been considered. However, 
explicit power series expressions are available |19|, although the coefficients in the series still involve 
determinants that need to be evaluated. 

Here we employ the approach of Ref. |24|, which is based on the relation to differential equations 
of Painleve type discovered by Jimbo at al. 1251 . to obtain series expansions or numerical estimates of 
these functions, compare also 1201 1271 1221 1231 for connections between random matrices and Painleve 
transcendents. In this paper, we generalize this approach to calculate series expansions for small spacings 
up to order 50 for the higher level-spacing distributions of the GUE. The expansions provide a detailed 
account of the small-spacing behavior of these distributions, and thus can be applied to analyze data. The 
exact coefficients involve ratios of large integer numbers, which makes them somewhat difficult to present 
in a paper. Therefore, the analytic expressions of the expansion coefficients are not included; they are 
provided in form of MATHEMATICA |26| code on the author's homepage |27|. The other two classical 
ensembles, the GOE and the GSE, can be treated in a similar fashion, although the calculation is somewhat 
more involved |24|. The results for these two universality classes will be published separately. 

Why are we interested in calculating these expansion to such high order in si There are several reasons. 
First of all, the data suggest that the series have infinite radius of convergence, thus the level-spacing 
distributions appear to be analytic functions of s. This means that, the more coefficients we know, the larger 
the interval where the truncated series gives a good approximation of the spacing functions. Secondly, the 
higher level spacing distributions have leading terms of increasing power in s, so in order to at least obtain 
the leading terms of the higher spacing function we have to expand all function during the calculation up 
to that given order in s. In our case, expanding up to order s 50 just suffices to gives us the leading terms 
of the spacing distribution of levels separated by five other levels. Thirdly, although it is, in principle, 
possible to integrate the system of coupled differential equations numerically and to obtain "numerically 
exact" results for the spacing functions, this does require the knowledge of at least the leading terms in the 
expansions as we cannot integrate directly from s = 0. Finally, it turns out to be very difficult to obtain a 
numerically exact distribution function in this way, because the distribution functions arise as differences 
between functions that grow exponentially with s, and thus requires high precision numerics. In fact, the 
author was not able to integrate numerically the set of equations for the spacing functions with four or five 
intermediate levels, even though the algebraic computer package MATHEMATICA 1 26 1 in principle allows 
for arbitrary precision arithmetics. 

The paper is organized as follows. After this introduction, the method used to calculate the spacing 
distribution is briefly summarized, following closely the discussion of Ref. 1241 . After that, the series 
expansion data are presented and their applicability is discussed. We end with some concluding remarks. 

2 Level-spacing distributions 

The basic objects that usually are considered are the probabilities E n (s) that an interval of length s contains 
exactly n eigenvalues. Here, s denotes the energy-level spacing in units of the mean level spacing. The gap 
probability Eq(s) that an interval of length s contains no eigenvalue at all can be expressed as a Fredholm 
determinant of a certain integral operator K with a sine kernel, see 1 10 1 and references therein. Jimbo et 
al. |25 1 showed that this Fredholm determinant can be expressed as 




(1) 
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where the function a(x; A), considered as a function of the variable x whereas A is regarded as a parameter, 
is a solution of a Painleve V differential equation 



( s ^0 +4 ( 



da 
x — a 

ax 



da 
x — a 

dx 



da\ 2 
dx ) . 



= 0. 



The parameter A is introduced via the boundary condition 

a(x; A) = x + o(x) 

7T 

for x — > 0. For the GUE, D(s; A) is the generating function of the probabilities E n (s), thus 
E n {s) = D n {s), 



where we defined 

D (s) 
D n (s) 



D(s;l) 

(-1)" d n D{s;X) 



n > 0. 



(2) 

(3) 

(4) 

(5) 
(6) 



The level-spacing distributions P n (s), which are the probability distributions to find two energy levels at a 
distance s with n levels lying in-between, are then obtained as 



Pn(s) 



ds 2 



{n - m + 1) E m (s) 



(7) 



m— 



where Pq(s) is the usual nearest-neighbour spacing distribution. 

Now, it is rather straightforward to derive series expansions for the distributions P n (s) by making a 
power-series ansatz for the solution a(x; A) of the differential equation 10 with the appropriate boundary 
condition (|3}- The distributions P n (s) were computed with the commercial algebraic computer program 
MATHEMATICA 1261 . Whereas it is rather easy to calculate the first few series coefficients in this way, a 
straightforward implementation of the equations in an algebraic programming language like MATHEMAT- 
ICA will generally not succeed to calculate higher orders. In order to achieve this, the calculation was split 
into smaller tasks and an iterative scheme was used, calculating the expansions coefficients one by one. In 
this way, it was possible to obtain the series expansions for the spacing functions P n (s) up to order s 50 . 



3 Series expansions 

Small-s expansions 



Pn(s) = ^Pn;kS k 



(8) 



of the level-spacing distributions P n (s) @ can be obtained as follows. As a first step, we define functions 

a n (x) by 



a (x) = <t(x;1), 

d n a(x; A) 



a n (x) 



9A* 



(9) 
(10) 
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The function cto (x) satisfies the differential equation (0 and behaves as 
1 



a (x) 



o{x) 



(11) 



for small x. This follows from Eq. (0, which also determines the small-a; behavior of a n (x) for n > 0, 
1 



<7l(x) 

a n {x) 



7T 

o{x) 



■ o(x), 
n > 1. 



(12) 
(13) 



For any N, the functions a n (x), < n < N satisfy a coupled set of differential equations that are 
obtained by taking the derivative of the original Painleve equation (0 with respect to A and putting A = 1. 
For instance, for N = 2 the set of differential equations becomes 



(xcjq) 2 + 4(xct - crn)(xa' - a + <t' 2 ), 
x 2 o'[o'q + 2(xa[ - ai)(xa' — a + a' Q 2 ) 
+2(xa' - <J )(xa[ - cr x + 2a[a' ), 
x 2 g'{gI + (xa'l) 2 

+2(xa' 2 — (T2)(x<t' - an + <r' a 2 ) 
+4(x<ri - o-i)(xfji - a% + 2a[a' ) 
+2(xa' Q — <TQ)(xa' 2 — (T 2 + 2a 2 (J o + 2cr^ 2 ). 



(14) 
(15) 



(16) 



In order to obtain series expansions for the level-spacing distributions P n {s) with < n < N, we first 
compute series expansions for the functions a n (x) for < n < N. This is done by inserting polynomials 



a n (x) 



K 

£< 

k=Q 



^ (9 



(17) 



of order K into the set of differential equations obtained from and implementing the initial conditions 
(fTTl-(tT3t by setting c„ ;0 = for < n < N, c 0; i = c\.\ = -1 and c„ : i = for 2 < n < N. The initial 
conditions ( II 11 and H2i motivate the convenient choice of the expansion variable x/ir rather than x. 

From the expansions of the functions a n (x), < n < N, we can derive expansions for D n (s) with 
< n < N by using Eqs. Q, ©, and (|6j. For example, we have 



exp 



D (s) 
Dx(s) = D (a) 



' Co (x) dx 



' o\ (x) dx 



(18) 



(19) 



D 2 (s) 



D (s) 



' (72 {x) dx 



+ 



' <Ji (x) dx 



(20) 



© 2003 WILEY- VCH Verkg GmbH & Co. KGaA. Weinhcim 



pss header will be provided by the publisher 



7 



Table 1 Numerical values of expansion coefficients ^8} of the level-spacing distributions Pq, Pi, P2, P3, Pi, and Pg. 
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and so forth. Finally, this translates into an expansion for the level-spacing distributions P n (s) via Eqs. (0} 
and 0. 

As Eq. Q involves two derivatives with respect to the variable s, we need to obtain the expansions 
of the functions a n (x) to two additional orders in x; so in order to achieve a result for P n (s) that is 
correct to order 50 in s we have to compute the expansions of a n (x) to order 52 in x/ir. Furthermore, we 
need to expand the exponential function in Dq(s) d 1 St t, whereas the integration on x is trivially performed 
on polynomials. As a result, the expansion coefficients p n -k (|U are expressions involving rather lengthy 
rational numbers and powers of it. The expansions for P n (s) with < n < 5 are given in AppendixlAl for 
higher n > 5 the lowest order in the small-s expansion of P n (s) is larger than 50. The numerical values of 
the expansion coefficients are given in Tabled compare also the numerical values of the coefficients po:fc 
for k < 32 given in Ref. |7 Table 4.1]. The modulus of the coefficients p n - t k appears to decrease relatively 
rapidly with increasing order k, so we may expect the series to converge on a rather large domain around 
the origin. 



4 Numerical integration 

To give an account of the accuracy of the expansions, we compare them to the result of a numerical 
integration of the coupled set of differential equations ( I14> -( I16> derived from the original Painleve equation 
0. These are also obtained using MATHEMATICA 1261 . In practice, we use initial values for the functions 
a n (x) for small, but non-zero x which are obtained from their expansions, and then integrate the system to 
larger values of x. 

In fact, it turns out that this is not as easily done as stated. The reason is that some of the functions 
involved become very small or very large as s is increased, and the computation of P n (s) essentially 
involves the cancellation of these large terms. In particular, the functions (J n (x) behave asymptotically like 

El 

a (x) ~ - l -x 2 (21) 
. , n\ exp(nx) 




Fig. 1 The spacing distribution Po(s) as obtained by numerical integration (solid line), the corresponding 
small-s expansion J23I and the asymptotic form 1301 . 
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Fig. 2 The spacing distribution Pi (s) as obtained by numerical integration (solid line), the corresponding 
small-s expansion J24I and the asymptotic form <31> . 



for large argument x. So, for n > 0, <r n (x) grows exponentially with x, whereas P n (s) is at most of 
order one, and eventually decreases as cxp(— 7r 2 s 2 /8) for large s, see Appendix 151 where the asymptotic 
behavior of the functions P n (s) for large spacings s is discussed. 

Results for n = 0, n = 1 and n = 2 are shown in Figs. ^ I3amd|3j respectively. Apparently, the 
numerical solutions agree well with the series expansions of AppendixlAlup to about s < 2.5. For n = 0, 
the distribution function is well reproduced by the asymptotic form of Eg. J30l > down to s > 1.5. Similarly, 
as shown in Figs.[5]and[3] the small-s expansions of Eg s . I24landl25 Idescribe the functions Pi (s) and P^ (s) 
for s < 2.5. The asymptotic forms of Eqs. <I3 1 i — d33l reproduce the functions Pi(s) and P2(s) very well 
for s > 3 and s > 4.5, respectively. In principle, it should be possible to extend the s _1 expansion of 
Eq. i31\ as well, at least for fixed values of n, and thus improve the situation considerably. 



expansion 
asymptotic 
numerical 




Fig. 3 The spacing distribution P2 (s) as obtained by numerical integration (solid line), the corresponding 
small-s expansion J25> and the asymptotic form <31> . 
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5 Concluding remarks 

We presented series expansions for the spacing distributions P n (s) of the GUE. The expansion was cal- 
culated up to order 50 in the spacing s, and for all distributions P n (s). In practice, this means that for 
< n < 5 the leading terms up to order s 50 were obtained, whereas for n > 5 the leading order is larger. 
The expansions agree with the complete series solutions obtained by Mehta 1 19 1, which were derived in a 
different way. 

The explicit expansions give a very precise account of the small-spacing part of the spacing distributions, 
and might prove useful for comparisons with experimental or numerical data. The leading asymptotic 
behavior for large spacing was obtained in Ref. [24 1 by similar means. 

Within the same framework, albeit slightly more involved, one can derive analogous series expansions 
for the spacing distributions of the GOE and the GSE. 

It is conceivable that other forms of expansions of the level spacing distributions P n (s) for small s might 
improve the convergence. A candidate might be to start from a Wigner surmise type form as discussed in 
Ref. ED- 
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A Leading terms of series expansions 

Here we list the leading terms of the small-s expansions for P„(s). The complete expansions can be 
downloaded as Mathematica code from |27|. The numerical values of the coefficients up to order s 50 are 
given in Tabled 

7T 2 S 2 2TT 4 S 4 7T 6 S 6 7T 6 S 7 2 TT 8 S 8 11tT 8 S 9 2 TT 10 S 10 

P (s) = 1 1 1 ... (23) 

w 3 45 315 4050 14175 496125 467775 

~ 4050 496125 + 13395375 170188239375 2679075000 + " ' 
P() tt 12 s 14 17tt 14 s 16 15777t 16 5 18 

2< ^ ~ 5358150000 1181472075000 + 2859162421500000 "' 

TT 20 S 23 13 7T 22 s 25 

3 ^ ~ 9378525331350000000 1702202347640025000000 + " ' ( 

_30 34 

PJs) = ... (27) 

4V ; 30645402510264863844600000000000 

TT 42 S 47 

Pais) = ... (28) 

by ' 255963589608666174754500410100972300000000000000 



For the level-spacing distributions P n (s) with n > 5 the leading terms are of higher order than s 



51) 



B Asymptotic behavior 

The large-s asymptotics of the level-spacing distributions P n (s) are also known. For n = 0, Dyson's 
asymptotic result 1 29 , 7 1 reads 
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and thus 




(30) 




(31) 



where 



B n — 2 




TT 



(32) 



rn—1 



From this, we can easily calculate the coiTesponding asymptotic behavior P^ a) (s) as 



*•<•> = £ ± < 

m=0 



(33) 



which is just Eq. Q applied to the asymptotic expressions. The resulting terms are rather involved; how- 
ever, it is apparent from Eqs. ( I29> . (13 1 1 . and J33i that the leading behavior in all cases is exp(— 7r 2 s 2 /8). 
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